Nominal SARAL Transfer Function

This paper gives a calculation of the range correction and cross section of the SARAL (Satellite with Argos and ALtiKa) Indian/French ocean radar satellite retroreflector array assuming the cube corners are coated and have a dihedral angle offset of about 1.5 arcseconds to account for velocity aberration. The cubes are assumed to all have the same orientation within the mounting. The derived range correction may be applied in precise orbit determination analyses that use Satellite Laser Ranging (SLR) data to SARAL

Long-term variability of CO2 and O in the Mars upper atmosphere from MRO radio science data

We estimate the annual variability of CO2 and O partial density using approximately 6years of Mars Reconnaissance Orbiter (MRO) radio science data from August 2006 to January 2012, which cover three full Martian years (from the northern hemisphere summer of 28 to the northern hemisphere summer of 31). These two elements are the dominant species at the MRO periapsis altitude, constituting about 70-80% of the total density. We report the recovered annual cycle of CO2 and the annual and seasonal cycle of O in the upper atmosphere. Although no other observations are available at those altitudes, our results are in good agreement with the density measurements of the Mars Express Spectroscopy for Investigation of Characteristics of the Atmosphere of Mars, which uses stellar occultations between 60 and 130km to determine the CO2 variability, and with the Mars Global Reference Atmospheric Model 2010 for the O annual and seasonal variabilities. Furthermore, the updated model provides more reasonable MRO drag coefficients (CD), which are estimated to absorb mismodeling in the atmospheric density prediction. The higher content of dust in the atmosphere due to dust storms increases the density, so the CDs should compensate for this effect. The correlation between the drag coefficient and the dust optical depth, measured by the Mars Odyssey Thermal Emission Imaging System (THEMIS) instrument, increases from 0.4 to 0.8 with the a priori and adjusted models, respectively. The trend of CDs not only confirms a substantial improvement in the prediction of the atmospheric density with the updated model but also provides useful information for local dust storms, near MRO periapsis, that cannot be measured by the opacity level since THEMIS does not always sample the southern hemisphere evenly

VGOS VLBI Intensives between MACGO12M and WETTZ13S for the rapid determination of UT1-UTC

In this work, we present a status update and results of the designated research and development VLBI Intensive program VGOS-INT-S, observed between MACGO12M and WETTZ13S for the rapid determination of the Earth's phase of rotation, expressed via UT1-UTC. The main novelty of these sessions is the use of a special observation strategy, rapidly alternating between high- and low-elevation scans, enabling an improved determination of delays caused by the neutral atmosphere. Since 2021, 25 Intensive sessions have been observed successfully. In early 2022, VGOS-INT-S was among the most accurate Intensive programs with an average formal error $\sigma_{UT1-UTC}$ of 3.1 $\mu$s and a bias w.r.t. IERS C04 of 1.1 $\mu$s. Later, the session performance decreased due to multiple technical difficulties.Comment: 8 pages, 5 figure

Topography of the Moon from the Clementine Lidar

Range measurements from the lidar instrument carried aboard the Clementine spacecraft have been used to produce an accurate global topographic model of the Moon. This paper discusses the function of the lidar; the acquisition, processing, and filtering of observations to produce a global topographic model; and the determination of parameters that define the fundamental shape of the Moon. Our topographic model: a 72nd degree and order spherical harmonic expansion of lunar radii, is designated Goddard Lunar Topography Model 2 (GLTM 2). This topographic field has an absolute vertical accuracy of approximately 100 m and a spatial resolution of 2.5 deg. The field shows that the Moon can be described as a sphere with maximum positive and negative deviations of approx. 8 km, both occurring on the farside, in the areas of the Korolev and South Pole-Aitken (S.P.-Aitken) basins. The amplitude spectrum of the topography shows more power at longer wavelengths as compared to previous models, owing to more complete sampling of the surface, particularly the farside. A comparison of elevations derived from the Clementine lidar to control point elevations from the Apollo laser altimeters indicates that measured relative topographic heights generally agree to within approx. 200 in over the maria. While the major axis of the lunar gravity field is aligned in the Earth-Moon direction, the major axis of topography is displaced from this line by approximately 10 deg to the cast and intersects the farside 24 deg north of the equator. The magnitude of impact basin topography is greater than the lunar flattening (approx. 2 km) and equatorial ellipticity (approx. 800 m), which imposes a significant challenge to interpreting the lunar figure. The floors of mare basins are shown to lie close to an equipotential surface, while the floors of unflooded large basins, except for S.P.-Aitken, lie above this equipotential. The radii of basin floors are thus consistent with a hydrostatic mechanism for the absence of significant farside maria except for S.P.-Aitken, whose depth and lack of mare require significant internal compositional and/or thermal heterogeneity. A macroscale surface roughness map shows that roughness at length scales of 10(exp 1) - 10(exp 2) km correlates with elevation and surface age

Towards the GEOSAT Follow-On Precise Orbit Determination Goals of High Accuracy and Near-Real-Time Processing

The US Navy's GEOSAT Follow-On spacecraft (GFO) primary mission objective is to map the oceans using a radar altimeter. Satellite laser ranging data, especially in combination with altimeter crossover data, offer the only means of determining high-quality precise orbits. Two tuned gravity models, PGS7727 and PGS7777b, were created at NASA GSFC for GFO that reduce the predicted radial orbit through degree 70 to 13.7 and 10.0 mm. A macromodel was developed to model the nonconservative forces and the SLR spacecraft measurement offset was adjusted to remove a mean bias. Using these improved models, satellite-ranging data, altimeter crossover data, and Doppler data are used to compute both daily medium precision orbits with a latency of less than 24 hours. Final precise orbits are also computed using these tracking data and exported with a latency of three to four weeks to NOAA for use on the GFO Geophysical Data Records (GDR s). The estimated orbit precision of the daily orbits is between 10 and 20 cm, whereas the precise orbits have a precision of 5 cm

The Precise Orbit and the Challenge of Long Term Stability

The computation of a precise orbit reference is a fundamental component of the altimetric measurement. Since the dawn of the modern altimeter age, orbit accuracy has been determined by the quality of the GPS, SLR, and DORIS tracking systems, the fidelity of the measurement and force models, and the choice of parameterization for the orbit solutions, and whether a dynamic or a reduced-dynamic strategy is used to calculate the orbits. At the start of the TOPEX mission, the inaccuracies in the modeling of static gravity, dynamic ocean tides, and the nonconservative forces dominated the orbit error budget. Much of the error due to dynamic mismodeling can be compensated by reduced-dynamic tracking techniques depending on the measurement system strength. In the last decade, the launch of the GRACE mission has eliminated the static gravity field as a concern, and the background force models and the terrestrial reference frame have been systematically refined. GPS systems have realized many improvements, including better modeling of the forces on the GPS spacecraft, large increases in the ground tracking network, and improved modeling of the GPS measurements. DORIS systems have achieved improvements through the use of new antennae, more stable monumentation, and of satellite receivers that can track multiple beacons, and as well as through improved modeling of the nonconservative forces. Many of these improvements have been applied in the new reprocessed time series of orbits produced for the ERS satellites, Envisat, TOPEX/Poseidon and the Jason satellites, and as well as for the most recent Cryosat-2 and HY2A. We now face the challenge of maintaining a stable orbit reference for these altimetric satellites. Changes in the time-variable gravity field of the Earth and how these are modelled have been shown to affect the orbit evolution, and the calibration of the altimetric data with tide gauges. The accuracy of the reference frame realizations, and their projection into the future remains a source of error. Other sources of omission error include the geocenter for which no consensus model is as of yet applied. Although progress has been made in nonconservative force modeling through the use of detailed satellite-specific models, radiation pressure modeling, and atmospheric density modeling remain a potential source of orbit error. The longer term influence of variations in the solar and terrestrial radiation fields over annual and solar cycles remains principally untested. Also the long term variation in optical and thermal properties of the space vehicle surfaces would contribute to biases in the orbital frame if ignored. We review the status of altimetric precision orbit determination as exemplified by the recent computations undertaken by the different analysis centers for ERS, Envisat, TOPEX/Poseidon, Jason, Cryosat2 and HY2A, and we provide a perspective on the challenges for future missions such as the Jason-3, SENTINEL-3 and SWOT

SLR Station Recovery, Center of Frame Motion, and Time Varying Gravity

Weekly station position estimates, beginning with 1993, are derived from the ITRF2008-based SLR processing of up to four satellites: Lageos 1, Lageos2, Starlette, and Stella. Helmert parameters obtained from c omparison of weekly SLR station positions and the a-priori SLRF2008 station complement are evaluated for geocenter motion and scale. Two me thods for modeling time varying gravity are employed in the SLR satel lite POD processing, with GGM03S serving as the static gravity field. Both methods forward model atmosphere gravity derived from 6-hour ECM WF pressure data. The standard approach applies an annual 20x20 field estimated from 4 years of GRACE data, and the IERS2003 recommended linear rates for C20, C30, C40, C21, and S21. The alternate approach us es a new set of low-order/degree 4x4 coefficients estimated weekly fr om SLR & DORIS processing to 10 satellites from 1993-2012. This exper imental tvg4x4 model has been shown to improve the TOPEX, Jason-1, and Jason-2 altimeter satellite orbits,. In this paper we apply the more detailed time-variable gravity modeling to the SLR satellite POD pro cessing and subsequent reference frame analyses. For this study we will evaluate the orbit differences (periodic and secular) for the satel lites concerned, characterize the impact on the station coordinate solutions, and the impact on reference frame parameters (geocenter and s cale)

Time Variable Gravity modeling for Precise Orbits Across the TOPEX/Poseidon, Jason-l and Jason-2 Missions

Modeling of the Time Variable Gravity (TVG) is believed to constitute one of the the largest remaining source of orbit error for altimeter satellite POD. The GSFC operational TVG model consists of forward modeling the atmospheric gravity using ECMWF 6-hour pressure data, a GRACE derived 20x20 annual field to account for changes in the hydrology and ocean water mass, and linear rates for C20, C30, C40, based on 17 years of SLR data analysis (IERS 2003) using the EIGEN-GL04S1 (a GRACE+Lageos-based geopotential solution). Although the GSFC Operational model can be applied from 1987, there may be long-term variations not captured by these linear models, and more importantly the linear models may not be consistent with more recent surface mass trends due to global climate change, We have evaluated the impact of TVG in two different wavs: (1) by using the more recent EIGEN-6S gravity model developed by the GFZ/GRGS tearm, which consists of annual, semi-annual and secular changes in the coefficients to 50x50 determined over 8(?) years of GRACE+Lageos+GOCE data (2003-200?): (2) Application of 4x4 solutions developed from a multi satellite SLR+DORIS solution based on GGM03S that span the period from 1993 to 2011. We have evaluated the recently released EIGEN6s static and time-varying gravity field for Jason-2 (J2). Jason-I (J1), and TOPEX/Posiedon (TP) Precise Orbit Determination (POD) spanning 1993-2011. Although EIGEN6s shows significant improvement for J2POD spanning 2008 - 2011, it also shows significant degradation for TP POD from 1992. The GSFC 4x4 time SLR+DORIS-based series spans 1993 to mid 2011, and shows promise for POD. We evaluate the performance of the different TVG models based on analysis of tracking data residuals use of independent data such as altimeter crossovers, and through analysis of differences with internally-generated and externally generated orbits

Ice mass change in Greenland and Antarctica between 1993 and 2013 from satellite gravity measurements

We construct long-term time series of Greenland and Antarctic ice sheet mass change from satellite gravity measurements. A statistical reconstruction approach is developed based on a Principal Component Analysis to combine high-resolution spatial modes from the Gravity Recovery and Climate Experiment (GRACE) mission with the gravity information from conventional satellite track-ing data. Uncertainties of this reconstruction are rigorously assessed; they include temporal limitations for short GRACE measurements, spatial limitations for the low-resolution conventional tracking data measurements, and limitations of the estimated statistical relationships between low and high degree potential coe�cients re ected in the PCA modes. Trends of mass variations in Greenland and Antarctica are assessed against a number of previous studies. The resulting time series for Greenland show a higher rate of mass loss than other methods before 2000, while the Antarctic ice sheet appears heavily in uenced by interannual variations

Solar System Expansion and Strong Equivalence Principle as Seen by the NASA MESSENGER Mission

The NASA MESSENGER mission explored the innermost planet of the solar system and obtained a rich dataset of range measurements for the determination of Mercury's ephemeris. Here we use these precise data collected over seven years to estimate parameters related to General Relativity and the evolution of the Sun. These results confirm the validity of the Strong Equivalence Principle with a significantly refined uncertainty of the Nordtvedt parameter eta=(-6.6 plus or minus 7.2)x10(exp -5) By assuming a metric theory of gravitation, we retrieved the Post-Newtonian parameter beta = 1 + (-1.6 plus or minus 1.8)x10(exp -5) and the Sun's gravitational oblateness, J(sub 2 solar)=(2.246 plus or minus 0.022)x10(exp -7). Finally, we obtain an estimate of the time variation of the Sun gravitational parameter, G (raised dot)solar mass/G solar mass =(-6.13 plus or minus 1.47)x10(exp -14), which is consistent with the expected solar mass loss due to the solar wind and interior processes. This measurement allows us to constrain |G(raised dot)|/G to be less than 4 x 10(exp -14) yr(exp -1)